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Abstract. A nonlinear model representing the tribological problem of a thin solid lubricant 
layer between two sliding periodic surfaces is used to analyze the phenomenon of hysteresis at 
pinning/depinning around a moving state rather than around a statically pinned state. The 
cycling of an external driving force F ex t is used as a simple means to destroy and then to 
recover the dynamically pinned state previously discovered for the lubricant center-of-mass 
velocity. De-pinning to a quasi-freely sliding state occurs either directly, with a single jump, or 
through a sequence of discontinuous transitions. The intermediate sliding steps are reminiscent 
of phase-locked states and stick-slip motion in static friction, and can be interpreted in terms 
of the appearance of travelling density defects in an otherwise regular arrangement of kinks. 
Re-pinning occurs more smoothly, through the successive disappearance of different travelling 
defects. The resulting bistability and multistability regions may also be explored by varying 
mechanical parameters other than -F oxt , e.g. the sliding velocity or the corrugation amplitude of 
the sliders. 



1. Introduction 

Nonlinear systems driven far from equilibrium exhibit a very rich variety of complex spatial and 
temporal behaviors pQ. In particular, in the emerging field of nanoscale science and technology, 
understanding the nonequilibrium dynamics of systems with many degrees of freedom which 
are pinned in some periodic potential, as is commonly the case in solid-state physics, is often 
becoming an issue. Friction belongs to this category too, because the microscopic corrugation of 
the mating surfaces may interlock j2j[3]. Simple phenomenological models are important, as they 
often give not only qualitative understanding of experimental findings, but also fair quantitative 
agreement with nanoscale tribology data, and with realistic simulations of sliding phenomena [I]. 
In this line of simplified approaches, studies are typically restricted to describing microscopic 
dynamics in one (ID) or two (2D) spatial dimensions. The substrates defining the moving 
interface are modelled in a simplified way as purely rigid surfaces or as one- or two-dimensional 
arrays of particles interacting through simple (e.g., harmonic) potentials. Despite such a crude 
level of description, this class of approaches frequently reveals the ability of modelling the main 
features of the complex microscopic dynamics, ranging from regular to chaotic motion [5l[6j[7]. 

One of the pervasive concepts of modern tribology - with a wide area of relevant practical 
applications as well as fundamental theoretical issues - is the idea of free sliding connected 



with incommensurability. When two crystalline workpieces with incommensurate or misaligned 
lattices are brought into contact, the minimal force required to achieve sliding, i.e. the static 
friction, should vanish, at least provided the two crystals are stiff enough. In such a geometrical 
configuration, the lattice mismatch can prevent interlocking of the two periodic corrugations and 
the resulting collective stick-slip motion of the interface atoms, with a consequent dramatically 
reduced frictional force. Experimental observation of this sort of superlubric and anisotropic 
regime of motion has been reported recently [8, 9j . The paradigm of frictionless sliding is realized 
naturally by the ID Frenkel-Kontorova (FK) model (see Ref. [S] and references therein). However 
the physical contact between two solids is generally mediated by so-called "third bodies", and 
the role of incommensurability has been recently extended [TU] in the framework of a driven ID 
model inspired by the tribological problem of two sliding interfaces with a thin solid lubricant 
layer in between. The frictional interface is thus characterized by three inherent length scales 
along the sliding direction: the periods of the bottom and top substrates, and the period of 
the embedded solid lubricant structure. In particular, in the presence of a uniform external 
driving velocity, the interplay of these incommensurate length scales can give rise to intriguing 
dynamical phase locking phenomena and surprising velocity quantization effects [12] . 

Previous numerical and theoretical studies of this confined tribological model |11[ [T2l [13] 
discovered a quantization of the lubricant center-of-mass (CM) relative velocity and found it to 
be related to the pinning of topological density excitations (kinks) to the substrate of closest 
periodicity. More recent work [JJ] highlighted a strict analogy of these dynamical pinning 
phenomena to the ordinary commensurate pinning of static friction [15} 116] . The proposed 
mapping between this dynamical pinning and that of static friction was explored numerically by 
analyzing the effect of an additional external driving force F e ^ t , equal for all lubricant particles. 
Dynamical pinning is signified by the lubricant CM relative velocity remaining robustly locked to 
the quantized plateau value (a value strictly and analytically determined by spatial periodicity 
ratios alone) up to a critical force threshold, above which quantization is destroyed. 

It was also found that as long as inertial effects are non negligible compared to dissipative 
forces (underdamped regime of motion), the adiabatic variation (increase and decrease) of the 
external driving force gives rise to a large hysteresis loop in the v cm - F ex t characteristics, not 
unlike depinning in static friction [5[ [15] . The present paper focuses precisely on the hysteretic 
behavior around a dynamical quantized steady state that this system exhibits, and discuss 
similarities and differences between such a dynamical locking and the more usual static pinning. 
By exploiting configurations where the dynamics of individual kinks is easy to monitor visually, 
the mechanism of hysteresis will be clarified. Given the practical difficulty of an experimental 
setup where an equal driving force is applied to each lubricant particle on the fly, the F ext 
term is may be seen more as a useful mathematical device rather than a realistic suggestion 
for future measurements aimed at studying dynamical depinning. On the other hand, we will 
bring concrete examples of the hysteretical destruction and recovery of the CM velocity plateau 
by means of parameters other than F ex t being cycled. The cycling of the substrate sliding 
velocity or of the applied load sketch practical possibilities to address the dynamical hysteresis 
in experimental tribological investigations. 

2. Confined lubricant model: numerical simulations 

We will work with the one-dimensional generalization of the standard FK model introduced in 
Ref. [HIE], consisting of two rigid sinusoidal substrates, of spatial periodicity a + and a_, and a 
chain of harmonically interacting particles, of equilibrium length ao, mimicking the sandwiched 
lubricant layer, as sketched in the inset of Fig. [Tfl The two substrates move at a constant 

The harmonicity of interactions within the lubricant chain is merely a simplifying assumption, since test 
simulations with anharmonic inter-particle potentials (e.g. Morse and Lennard- Jones) also reveal the ubiquity 
of the observed phenomenology. 




Figure 1. (Color online) Normalized velocity of the center of mass, f C m/^ext) as a function 
of the chain stiffness KF + /a + , for t> ext = 0.1 (F+a+fm) 1 / 2 , 7 = 0.1 (i ? +m/a+) 1//2 , and 
r + = (1 + 7r 1//4 )/2. Crosses: one-to-one kink coverage 9 = 1 (r_ — 7.036); diamonds: 
commensurate kink coverage 8 = 5/4 (r_ ~ 8.795); circles: incommensurate kink coverage 
9 = (1 + 10 1//2 )/3 (r_ ~ 9.762). Dashed line: the quantized-plateau velocity ratio of Eq. 
Note the logarithmic scale in the abscissa. Inset: a sketch of the driven 3-lengthscale confined 
model. 



relative velocity v cx t = V- — v+. In particular, we select the reference frame where v+ = and 
u_ = f ext . The equation of motion of the i-th lubricant particle is: 



1 

mxj = — 
2 



F + sm — Xi + F_ sm — (Xj — v ext t) 



+K(x i+ i + Xi-\ - 2xi) - 2^{xi - Vyt) + F ext , (1) 

where m is its mass. F± are the amplitudes of the forces due to the sinusoidal corrugation of the 
substrates. By default, we set F_/F+ = 1 as the least biased choice, but we will explore the effect 
of modifying F_ in Sect. [3] below. K is the chain spring constant defining the harmonic nearest- 
neighbor interparticle interaction. The penultimate damping term in Eq. ([1]) originates from two 
symmetric frictional contributions adding as follows: —7 (xj— v + )— 7 (xj — u_) = —27 (xj — ^v cxt ), 
where 7 is a viscous friction coefficient accounting phenomenologically for degrees of freedom 
inherent in the real physical system (such as substrate phonons, electronic excitations, etc.) 
which are not explicitly included in the model; this fixes the reference speed of the the dissipative 
term: v w = |u ex t El- In order to probe the strength of quantization, and eventually address 
hysteresis, an additional constant force F ext is applied to all chain particles and varied up and 
down adiabatically. The infinite chain size is managed - in the general incommensurate case 
- by means of periodic boundary conditions (PBC) and finite-size scaling [13] . We set overall 
a_|_ = 1, m = 1, and F + = 1 as basic dimensionless units, and express implicitly all mechanical 
quantities in terms of natural model units obtained as combinations of these three basic units 



2 We choose this value of the velocity u w = ^oxt to which dissipation refers as the least biased option. Different 
choices, equivalent to choosing different 7+ and 7- dissipation coefficients to the two substrates, would at most 
change the quantitative details of the velocity-plateau boundaries, but not the qualitative nature of results. 



As previously found [TTJ [T21 [T3J [H] the detailed behavior of the driven system in Eq. (TjQ) 
depends crucially on the relative (in) commensur ability of the substrates and the chain. The 
relevant length ratios are defined by r± = a±/ao; we assume r_ > min(r + , r^ 1 ), whereby the 
(+) substrate has the closest periodicity to the lubricant, the (— ) slider the furthest. Under 
rather general dynamical conditions, the lubricant slides with a quantized mean velocity v plateau 
relative to the (+) substrate. The plateau phenomenon was explained by the static pinning of 
the topological solitons (kinks) that the embedded chain forms with the (+) substrate, to the 
(— ) slider I13j. Specifically, the quantized-plateau lubricant velocity ratio 

^cm ^plateau 1 

Vcxt v ext r + ' 

is strictly a function of the lubricant coverage r + of the (+) substrate [TT], i.e. of the absolute 
density (r+ — of kinks. For antikinks, r + < 1, this density is negative, and so is t> p iateau _ 

namely the lubricant slides backwards Although the quantized-plateau velocity depends 

uniquely on r+, the plateau dynamical stability and extension depend crucially on the kink 
coverage 

e = aJ-±^± = r(l-±-\ (3) 

a+ V r+J 

of the (— ) substrate (for antikinks, 9 < 0). Concretely, as a function, e.g. of the spring stiffness 
K, the quantized plateau is very prominent in a range of K of the order unity but weakens and 
eventually terminates for stiffer chains (larger K values), Fig. [TJ The plateau destabilization 
is complete for a general irrational 9, while, under suitable conditions detailed below, the 
plateau can survive up to indefinitely large K for commensurate kink coverage (rational 9). 
The quantized velocity plateau is finally particularly robust for perfect one-to-one matching of 
the soliton and the (-) slider periodicities, 9 = 1 [Hj. To illustrate these three typical cases, 
we consider r + = (1 + 7r 1 / 4 )/2 ~ 1.166 and the three values r_ = 9 (1 — r^ 1 ) -1 ~ 7.036, 8.795, 
and 9.762, corresponding to the values 6 = 1, 6 = 5/4 = 1.250, and = (1 + W 1/2 )/3 ~ 1.387 
respectively. The choice of r + near unity, is especially advantageous compared to values like the 
golden mean (1 + \/5)/2 ~ 1.618 often used, because it gives rise to well-separated individual 
kinks, which allow a more transparent analysis of the dynamics. Many qualitative features 
discussed for the specific ratios r± considered here are in fact also found for general values of 
r±, and thus this specific choice of length ratios should not be considered especially restrictive, 
as long as a correct distinction of different commensuration property of 9, Eq. ([3]), is made. 

The equations of motion (TjQ) are integrated using a standard fourth-order Runge-Kutta 
algorithm. The system is initialized with the chain particles placed at rest at uniform separation 
do, and the top substrate is made slide at the imposed constant velocity V- = w cx t- For F ext = 
and a wide range of model parameters, after an initial transient the system reaches a steady state, 
where all dynamical quantities other than particle positions fluctuate but show no systematic 
drift. For wide ranges of parameters, exemplified in Fig.[T]by the spring stiffness K, the lubricant 
reaches the expected plateau state of normalized time-averaged velocity fpiateau/^ext — 0.142, 
Eq. ([2]), the same for the three geometries introduced above. 

Adiabatic upward and downward variation of the external force -F ex t is realized by changing 
F ext in small steps and letting the system evolve at each step for a time long-enough for all 
transient stresses to relax. This allows us to gauge the robustness of the plateau state as a 
function of the system parameters, e.g. of K. In order to determine the critical values of F^tt, 
where the plateau is abandoned and retrieved, and in particular to do that with great accuracy 
and a reasonably small number of separate simulations, we first increment F ex t in steps of 
0.01 F+, and then reduce the step width using a bisection scheme around the critical force. 



(2) 




Figure 2. (Color online) Hysteresis in the v cm — F ext characteristics for a confined chain of 
intermediate spring stiffness (K = 5F + /a+), and length ratios r+ = (1 +vr 1 / 4 )/2, 9 = 5/4 
(r_ ~ 8.795). The behavior is shown for fast (v cxt = 0.1, upper panel) and slow (t> ext = 0.01, 
lower panel) drive. Adiabatic increase and decrease of F ex t (in steps of 10~ 3 F + ) are denoted 
by crosses and circles, respectively. Characteristic hysteretic multi-step features appear. Here 
7 = 0.1(-F + m/a-|-) 1//2 , and a chain of N = 387 lubricant particles is simulated. 

3. Results 

For concreteness, we begin with the specific example 6 = 5/4, and pick an intermediate value of 
the chain stiffness K = 5 F + / a + , common to all plateaus of Fig. [TJ We start investigating the 
plateau destruction/recovery induced by varying the external force F ext through a sequence of 
adiabatic increases and decreases [I?] . The resulting CM velocity is displayed in Fig. [2] for two 
different external driving velocities v ex f A clear hysteretic loop emerges, with qualitatively 
similar features for high (upper panel) and low (lower panel) values of v ex t- Interestingly, 
and somewhat unexpectedly, the hysteretic regions are systematically broader for larger sliding 
velocities v ex t- We will return to this point later on. 

The exact plateau state implies a kind of dynamical incompressibility, namely identically 
null response to perturbations or fluctuations trying to deflect the CM velocity away from 
its quantized value. Indeed, as long as F ex t remains below a critical threshold it does 

perturb each individual single-particle motion, but has no effect whatsoever on v cm , which 
remains exactly pinned to the quantized value, as is indeed expected of an incompressible 
state. This behavior contrasts with all observed non-plateau sliding states, where v cm increases 
monotonically with -Fext- This plateau state is reminiscent of the pinned state of static friction, 
where a minimum force (the static friction force) is required to initiate the motion. Except that 
here in the starting "pinned" plateau state the lubricant chain particles are moving relative to 
both substrates. The sudden upward jump of v cm taking place at F ext = F+ * can thus be termed 
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Figure 3. (Color online) The (K,F ext ) phase diagram illustrating the unpinning-repinning 
transitions for r + = (l + 7r 1 / 4 )/2, 7 = 0.1, w ex t = 0.1, and for (a) one-to-one kink coverage 9 = 1; 
(b) commensurate kink coverage 9 = 5/4; (c) incommensurate kink coverage 9 = (1 + 10 1//2 )/3. 
The white areas have perfect plateau dynamics, the dotted region indicates quasi-free sliding. 
Simulations done with N = 387 particles for (a) and (b) and with N = 781 particles for (c). 

a dynamical depinning. The depinning transition line , appears as a "first-order" line, with 
a finite jump Av in the average v cm and a clear hysteretic behavior: as F ext is reduced back, the 
depinned state survives below down to a significantly smaller F^~^, where perfectly quantized 
plateau sliding is retrieved, as illustrated in Fig. [2l Several hysteretic loops are in fact observed 
in Fig. [2} a qualitatively similar multi-step behavior appears also for 9 = 1 and 9 = (l + W^ 2 )/3. 
We shall return below to the nature of these steps. The large- -F ex t quasi- free sliding regime is 
characterized by v cm increasing continuously, roughly proportionally to F ext /~f, and superposed 
to this general translational motion, by chaotic single-particle movements, contrasted to the 
periodic (9 = 1, 5/4) or quasi-periodic (9 = (1 + 10 1//2 )/3) individual-particle oscillations in the 
plateau state. 

The values of F^ , F+^, and Av are nontrivial functions of the parameters. Specifically, Fig. [3] 
reports the K dependency of these critical forces in the three considered cases. The values of the 
critical forces are remarkably similar for K < 4, while important differences are observed as the 
springs become stiffer. In particular, for unity coverage (9 = 1) the plateau is very stable and 
extends to very large K, as expected in a fully commensurate case, see Fig. E)(a). In contrast, 
for noninteger 9 the plateau becomes more fragile for large K. For commensurate 9 = 5/4 the 
plateau width decreases with some fast power law of K, and becomes numerically difficult to 
detect beyond K ~ 60. For incommensurate 9 = (1 + 10 1//2 )/3 instead, the plateau shrinks and 
disappears at finite K = ~ 24: no sign of a quantized plateau is detectable, e.g., for 



K = 25. This unequal behavior for commensurate/incommensurate coverage 9 is understood 
in terms of the mapping of the dynamical sliding model to the static FK model, which was 
established in Ref. [Hj. The hysteretic depinning transition is observed through a significantly 
wide if-range in all three cases, but the depinning mechanism differs in some important detail. 

3.1. Fully commensurate 9 = 1 

As illustrated in Fig. O for 9 = 1 the plateau extends to very large K, in a range of F ext 
of decreasing width oc if . if -1 describes precisely the asymptotic decrease of the sinusoidal 
interparticle distance modulation, residual after solitons overlap one another in the large if limit. 
For very large if, outside the right end of Fig. [3^a), the asymptotic values of this curve lie 
entirely in the negative- -F ex t domain. The explanation is that it takes a negative external force 
to compensate the positive average dissipative "wind" force F w = —2j(v crn — v w ) acting on each 
lubricant particle. On the plateau state this wind force amounts to 

2 - r + 

F w = -27 (plateau ~ %) = 7 ^cxt • (4) 

In the absence of the external driving -F ex t) the wind force alone is sufficient to disrupt the 
plateau at large if, where it is more fragile, as seen at the large- if side of Fig. [TJ However, once 
F w is compensated away, the 9 =1 quantized plateau extends to indefinitely large if. 

The next result concerns hysteresis, still at 9 =1. Depinning is discontinuous and hysteretical 
as exemplified in Fig. EJ but only up to a large but finite critical stiffness if = if* ~ 330. Near 
if* the bistability range F^ — F c + ^ closes up with a power law F^ — F^ = B (if* — if ) a , not 
unlike what was observed in previous work for the golden-mean ratio Ref. |14j . Above the critical 
stiffness, for if > if*, F+' = F^, the depinning transition is continuous and characterized by 
what appears to be a mean- field power law v cm — w p i a teau oc (F ext — F^) 1 / 2 . For if approaching 
if* from below, the plateau is abandoned through different mechanisms depending on the model 
parameters. In Ref. [16J it was found that re-pinning in the continuous sine-Gordon model 
proceeds first through a series of "cavity-mode" states, and then a series of kink-antikink wave 
train states, and a similar scenario is exhibited also by the discrete FK chain [15j . We find that 
analogous phenomena occur here for the repinning to the dynamical plateau, with defects in the 
kink lattice taking the place of the kink-antikink pairs of the single-chain FK model. 

For soft enough chains, individual kinks are visible and well distinct. For example Fig. 2] 
(decreasing F cxt ) illustrates the mechanism supporting deviations from the plateau for if = 5, 
the same value as Fig. [5J A kink vanishes at a (— ) lattice site and joins a second kink to form a 
mobile "bi-kink" . This extra density accumulation "binds" substantially less than a kink to the 
minima of the (— ) potential. The external force F ey ± acts on the bi-kink density lump and drags 
it along to the right. Contrary to the bi-kink, the site with a missing kink ("no- kink") remains 
pinned to the (— ) potential well, and is not dragged by the external force F ey ±. The moving 
bi-kink breaks the "quantized" motion by one single particle, and is responsible for displacing 
the lubricant CM velocity a little bit away from the exact v p i a teau- 

The number of bi-kink - no-kink pairs tends to increase rapidly with increasing F cxt — F^~^. 
The force F^ necessary to nucleate the first bi-kink - no-kink pair is sufficiently large to sustain 
an avalanche of more bi-kink - no-kink pairs after the first defect is nucleated. Trains of bi-kinks 
cross the chain, producing essentially chaotic motions of the single lubricant particles, provided 
that F cxt 3> F^K When, starting from this dislodged, or depinned state, F ext is gradually 
reduced, bi-kink - no-kink pairs annihilate, the number of these pairs reducing steadily with 
time. The discrete, integer nature of the defect-pair number originates the (gently sloping) 
discrete downward staircase steps in the hysteresis loop, generally similar to those shown in 
Fig. [2] (for a different 9). Since the discrete effect of the disappearance of a single defect-pair 
becomes negligible in the infinite-size limit, the observed multi-step structure appears to be 
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Figure 4. (Color online) Soft chain {K = 5). Four snapshots of a 60-particles section of the 
lubricant chain and (+/— ) substrates (lower/upper sinusoids), at successive times separated by 
14 time units (a+ m/ F + ) 1 / 2 . The horizontal direction represents distance, the dots the particle 
positions Xj. Vertical displacements of dots measure the distance Xi — Xi-\ of a given particle 
to its left neighbor: on this scale, the horizontal solid and dashed lines indicate the average 
interparticle distance ao, and the (+) lattice parameter a+ respectively. The snapshots refer to 
r+ = (1 + vr 1 / 4 )/2, 9 = 1, F ext = 0.08136 (decreasing), and illustrate the crossing of the critical 
line F c + ^, with the recovery of the plateau state (see Fig. [2]) occurring through the disappearance 
of the last bi-kink - no-kink defect. The other parameters are 7 = 0.1, and v ex t = 0.1. This 
annihilation of a bi-kink against a stationary no-kink is best illustrated by the online animation 
repinning_thetal_K5.gif, which spans 70 time units, starting 11 time units before frame (a) and 
ending 17 time units after frame (d) of the present figure. For this and all animations we select 
the reference frame where the (— ) substrate, and thus all pinned kinks, are stationary. 

merely a finite size artifact, and for all that we can tell at present the infinite system should 
exhibit no staircase steps. In the depinned state, so long as F ext is strong enough, a bi-kink 
encounters a no-kink, interacts briefly, and then continues to travel. When instead F cxt is 
reduced below F^*, as in Fig. U the encounter of a bi-kink and a no- kink leads to reciprocal 
annihilation. The amplitude oscillation still visible (but quickly damped) at the right end side 
of the last frame of Fig. [J] reflects the waves dissipating the excess ( "binding" ) energy of the 
bi-kink - no-kink pair, in the process of recovering the perfect kink lattice. When finally the 
kink lattice gets rid of the last defect pair, the perfect plateau state is re-gained. 

For a stiff enough chain, individual kinks become spatially broad, and will for a fixed density 
extend over a size larger than the average inter-kink distance a+/(r+ — 1). In this limit the kink 
lattice reduces to a weak sinusoidal deformation, of amplitude oc K~ l superposed to the average 
interparticle density. Despite this difference with the strong kink lattice of the soft-chain case, 
the external-force-induced departure from the quantized velocity plateau occurs here through 
a mechanism similar to that illustrated above for the soft-spring case. A chain slippage by one 
particle (i.e. a distance ao) is promoted by a bi-kink and a no-kink moving in opposite directions: 
when they collide, the bi-kink - no-kink pair takes the aspect of a broad locally flat region of 
denser-than-average and less-dense-than-average lubricant in the otherwise perfect pinned kink 
lattice. As illustrated in Fig. [5ja), a local flattening defect forms in the soliton lattice, similar to 
the local amplitude suppression of a dragged charge density wave (CDW) [THJ [19]. This defect 
is characterized by a smooth "charge" separation, with the denser region being driven to the 




Figure 5. (Color online) Stiff chain {K = 50). Four successive snapshots of the substrates and 
lubricant chain, separated by time intervals of 9 model units. All notations and parameters are 
the same as in Fig. UJ except for K = 50, -F ex t = 0.00685 (decreasing), which falls in the region 
immediately above the critical line F^~^, before the recovery of the plateau state. The complete 
collision of the right-traveling bi-kink and left-traveling no-kink is best illustrated by the online 
animation unpinned_thetal_K50.gif, which spans 50 time units, starting 16 time units before 
frame (a) and ending 7 time units after frame (d). 

right and the more rarefied region to the left by the driving force, the external force acting like 
an electric field on a CDW insulator. These defects travel in opposite directions, as expected of 
a opposite charges driven by an electric field. The crucial difference with the soft-spring case 
(where as shown by Fig. [H the no-kink defect remains pinned to the (— ) lattice) is that here 
both defects, the by-kink and the no-kink, are mobile and dragged by the external force. As 
the two defects move apart, a perfect soliton lattice re-forms in between, Fig. E|d). In time, a 
right-moving bi-kink encounters a left-moving no-kink: these defects may again cross, or else 
they may bind and annihilate in pairs. Annihilation occurs when F ex t is reduced below -F c + ^, as 
in the soft-chain case of Fig. SJ When instead, F cxt > F^ the pair separates again, with the 
rightward "positive" and leftward "negative" flattenings suffering some phase shift, but traveling 
on, as in Fig.[5j As soon as all defects annihilate, the kink lattice is perfect, and the CM velocity 
recovers f p i a teau exactly. If the defect pairs form at regular spatial separation within the chain 
(with periodic boundary conditions) the corresponding moving pattern leads to time-periodic 
fluctuations of the CM velocity; that can also be seen as type-I intermittencies [20] . Otherwise, 
when defect motion is chaotic, an irregular CM dynamics is observed. For indefinitely growing 
chain stiffness K, each defect pair flattening region grows in size, eventually covering the entire 
finite-size simulation, which becomes at that point a poor representation of the infinite-size 
thermodynamical limit. 

Figure [6] draws the plateau boundaries relative to -F ex t> for varied external driving w cxt , for 
a rather stiff chain {K = 50). As the friction-drag reference force F w grows linearly with f ex t, 

it /I 

and this introduces a trivial compensating trend F c oc — F w , it is convenient to remove the 
appropriate linear drift by adding F w , Eq. to the critical forces. The static limit v ext = 
is smooth, and this indicates a regime of continuity from the static quasiperiodic 3-lengthscale 
model of Ref. [21] to the dynamical sliding. Strikingly, the plateau robustness against the 
external perturbing force F ext and the widths of the hysteretical regions benefit of increased 
driving speed. For large v ext — 7 the plateau destabilizes suddenly and eventually disappears. 
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Figure 6. (Color online) Driving- velocity dependency of the dynamical depinning and repinning 
forces F+* F+^, (shifted upward by the trivial F w oc v cxt contribution, Eq. Q). v ext is 

measured in model units of (F+a+fm) 1 / 2 ; the chain is rather hard (K = 50); r + = (1 + 7r 1 / 4 )/2, 
9 = 1 (r_ ~ 7.036), and 7 = 0.1. 




Figure 7. (Color online) Typical plateau arrangements of the 8 = 5/4 commensurate soft 
K = 5 (a) and hard A" = 50 (b) chain: a regular arrangement of bi-kinks (one every four kinks). 
The conventions and all other parameters are the same as in Fig. 01 but for F ext = —F w . 

3.2. Commensurate 8 = 5/4 

Having explored at length the 8 = 1 commensurability, we now turn to another kink lattice / slider 
system, still commensurate but with 8 = 5/4, a weaker commensurability than 9 = 1. At 

8 = 5/4, in the perfect-plateau state one kink out of four turns into a bi-kink, as illustrated 
in Fig. [7J (The bi-kinks of the present 8 > 1 case would be replaced by no-kinks for 8 < 1). 
The pre-existence of a regular array of such defects of the kink lattice allows for a significantly 
different depinning mechanism, compared to the totally commensurate 8 = 1 case. Defects of 
the kink lattice are already present prior to turning on the external force F ext , which only sets 
them into motion, without a need to create them. For soft springs, Fig. [7J(a), where the pinning 
energy barrier of these defects is large, Fig. [3] shows that the critical forces needed to set the 
defects into motion in this 9 = 5/4 case are very similar to those for 8 = 1. For harder springs, 
defects increase in size and affect several neighboring kinks now, as illustrated in Fig. W[b). 
These extended disturbances possess a much smaller pinning energy to the (— ) potential. As a 
consequence, the plateau state is now exceedingly weak, confined to an extremely narrow force 
range around — F w , see Fig. [3j The ordered arrangement of defects still warrants some amount 
of pinning, but the width F^ — F~^ of the pinned region decreases much faster than in the 

9 = 1 case as soon as the defect size exceeds the typical inter-defect distance a + /(r + — l)/9, 
here occurring for K ~ 10. 
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Figure 8. (Color online) 9 = (1 + 10 1//2 )/3 incommensurate soft chain K = 5 (a) and hard 
chain K = 50 (b): irregular alternation of kinks and bi- kinks. Pinning is realized for the soft 
chain, while even with F ext = — F w , so that v crn ~ Vpiateam the K = 50 hard chain is unpinned, 
with the defects slowly drifting along. The conventions and all other parameters are the same as 
in Fig. SI The online file unpinned_incommensurate.gif provides an animation of the situation 
of snapshot (b). A second multimedia file, unpinnedjncommensurate_drifting.gif, shows the 
defects drifting under the effect of a deviation of F ext by 10~ 3 force units in excess of — F w . 



3.3. Incommensurate 9 = (1 + 10 1//2 )/3 

Finally, at irrational 9 = (1 + 10 1/2 )/3, some kinks are replaced by bi-kinks, but the 
incommensuracy of the coverage leads to their irregular arrangement, as illustrated in Fig. [8j 
For a sufficiently soft chain, (represented by K = 5 in Fig. E](a)), the irregular distribution 
of single kinks and bi-kinks remains statically pinned to the minima of the (— ) substrate, 
with a finite barrier to overcome for a bi-kink to migrate to the next minimum. This barrier 
guarantees the existence and robustness of the CM quantized velocity plateau (with a first- 
order hysteretical boundary) in the present incommensurate case, pretty much like for the 
commensurate cases. This energy barrier protects the plateau against the movement of bi-kinks 
until K < ^Aubry — ^4. In contrast, for a harder chain (K > -f^ubry)' illustrated by K = 50 
in Fig. E](b) , the irregular distribution of single kinks and bi-kinks drifts through the chain at a 
speed approximately proportional to F ex i-\-F w , with no sign of any pinned plateau: this indicates 
that the energy barrier is here entirely removed by the irregular bi-kink configuration produced 
by incommensuracy. The kink-kink repulsion makes the bi-kinks increasingly extended objects 
as K increases, until they become so broad that crossing the maxima of the (— ) potential costs 
negligible energy: the bi-kink in the central region of Fig. EKb) exemplifies precisely one such 
slow hopping process. The transition between the soft-chain dynamically pinned regime and 
the stiff-chain fully unpinned state is analogous to the Aubry transition observed in the static 
situation described by the FK model. The kinks of the dynamical model play the role of the 
particles of the static model. 



3.4- Hysteresis when cycling other parameters 

By analogy to the single-chain FK model, cycling the external force F ext is conceptually the most 
natural way to abandon and recover, often hysteretically, the dynamical plateau. However, 
in practice, the experimental realization of a uniform force acting equally on each lubricant 
particle in flight is not trivial. On the other hand, the plateau can be abandoned and recovered, 
even when different parameters are cycled. Within the present model, the reason is that the 
dissipation 7-term has itself the effect of diverting the CM velocity away from v p i a teau- In a 
concrete laboratory configuration moreover, beside dissipative effects, other interactions too will 
tend to push the lubricant slide at speeds other than i> p i a teau- As an example, defects and 
grain boundaries will tend to pin statically the lubricant to either substrate [22J. These other 
"external" forces compete with the tendency to dynamical pinning: the latter tuned by other 
parameters, namely, in the language of our model, K F + F- and v ex t- Thus in a practical 



F 

0.5 0.6 0.7 0.8 0.9 " 1 1.1 1.2 1.3 1.4 




o— o increase 

nl I i I i I i I i I i I 

0.026 0.028 0.03 0.032 0.034 0.036 
F 

Figure 9. (Color online) Hysteresis loops found as the corrugation of the (— ) substrate F- is 
cycled down from (crosses) and back up to (circles) its value F + used in all other calculations, 
(a) At K = 160, near the plateau edge of Fig. [JJ it takes a small decrease in F_ to leave the 
plateau, while (b) when the plateau is very robust {K = 5), nonhysteretic depinning is observed 
for a corrugation amplitude F- far below unity. Simulations for r + = (1 + 7r 1 / 4 )/2, 9 = 1 
(r_ ~ 7.036), 7 = 0.1, v ext = 0.1. 

straightforward experiment, cycling quantities such as the the sliding speed, or the load applied 
to the sliders should lead to leaving/recovering the plateau dynamics, with hysteretic cycles 
similar to those exemplified by Fig. [2j 

To illustrate this point within our model, Fig. [9] depicts a first example of such a hysteretic 
cycle, where the load applied to the sliders, proportional to the upper slider corrugation 
is cycled. The plateau is abandoned hysteretically when F_ is decreased below critical values 
which depend strongly on the robustness of the pinned state, which is, in turn, a function of K 
and other model parameters. 

Along a similar scheme, the perfectly legitimate interpretation of Fig. [6] as a phase diagram 
suggests that the first-order line separating the free-sliding regime from the perfect plateau could 
be crossed by cycling f ext rather than F cxt . This cycle corresponds to tracking up and down 
the .Fext = dashed path drawn in Fig. The resulting loop, shown in Fig. [TUT a), depicts 
the expected bistability: v cxt is cycled up and down, and the perfect plateau is abandoned at 
much larger speed than where it is recovered. At large speed F w increases, the dissipative term 
dominates and makes the lubricant speed approach v w . 

The depinning transition may also occur continuously, when the transition line is crossed 
beyond the tricritical point, i.e. for K > K*, in the strongly dissipative region, where the 
the viscous damping rate r y/m is much larger than the vibrational frequencies, decreasing 
proportionally to K , of the soft kink lattice around the minima of the (— ) potential. In this 
regime the dynamical depinning is apparently second order. In this overdamped regime, shown 
for example in Fig. [TOT b). the forward and backward trajectories become indistinguishable, and 
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Figure 10. (Color online) (a) Hysteresis loop at the plateau edge in the f cm — ^ext characteristics 
for a confined chain of length ratios r+ = (1 + vr 1//4 )/2, 8 = 1 (r_ = 7.036). Adiabatic increase 
and decrease of v cxt are denoted by crosses and circles, respectively. Here 7 = 0.1 (F + m/a + ) 1 / 2 
and F cxt = 0, which corresponds to the dashed path of Fig. [6j (b) No hysteresis is observed in 
the overdamped regime (7 = 1.0 (F+m/ a^) 1 / 2 ) along the same path. 



hysteresis disappears. In this strongly dissipative regime, we find, instead of the hysteretic 
jumps, a nonlinear dependency of v cm versus the model parameters (here v ex t, but cycling F ext , 
F_, or K would lead to perfectly analogous results), without any bistability phenomena. 

4. Discussion and conclusions 

We have shown that starting from the quantized sliding plateau state, previously found for a 
simple tribological model of a confined layer, the sliding dynamics of the lubricant layer exhibits 
a large hysteresis when an additional external driving force F ex t trying to push v cm away from its 
quantized value is cycled. In analogy to depinning in ordinary static friction [15], the hysteretic 
dynamical behavior depends strongly on whether the system degrees of freedom have sufficient 
inertia (underdamped regime) or if, on the contrary, inertia is negligible (overdamped regime). 
Hysteretic versus continuous depinning occurs depending on whether the unpinning transition is 
crossed below or above a tricritical point where hysteresis close, and which marks the separation 
between the underdamped and the overdamped dynamics. 

Hysteresis arises due to the great robustness of the quantized dynamics, setting a large critical 
threshold F^ to the formation of mobile defects (initially depinned bi- kinks or no- kinks). Once 
at least one of these defects forms, an avalanche process leads to a discontinuous jump to a free 
or quasi-free sliding regime. Starting from the unpinned states, the plateau recovers only at 
a much smaller threshold F c + ^, representing the minimum driving force needed to sustain the 
motion of pre-existing mobile defects. 



Nontrivial differences with static friction occur. The first is that the dynamical pinning 
hysteresis cycle may be larger in situations where pinning itself could be intuitively considered 
more fragile, e.g., for larger external velocity. Another feature (presently under investigation, 
not discussed above) is that the sudden application of an external force can sometimes leave v cm 
locked to the quantized value, even if the applied force is larger than the dynamic depinning 
threshold F+T obtained instead through the adiabatic procedure sketched above. Once again, 
this is different from static depinning, usually requiring smaller force (than the static friction 
F s ) if applied suddenly [5]. 

The present study concentrates on zero temperature. At finite temperature, the energy barrier 
to the formation of defects such as bi-kinks and for defects "hopping" to neighboring pinning 
sites can be traversed by means of random thermal excitations. This means that at sufficiently 
low temperature the dynamical pinning should not change much. Even the hysteresis should 
remain, provided that parameters such as -F ex t are cycled much faster than the characteristic 
thermal relaxation times. Thermal effects are currently under closer investigation. 
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